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Abstract. The origin of the local population of young, cooling neutron stars is investigated with a population 
synthesis model taking into account the contribution of neutron stars born in the Gould Belt, in addition to those 
originating in the Galactic disk. We estimate their emission in the soft X-ray band as a function of distance and 
age and construct the Log N - Log S distribution. It is shown that the inclusion of neutron stars from the Gould 
Belt provides a good fit to the observed Log N - Log S distribution. As the Sun is situated inside the Gould 
Belt, one can naturally explain the comparative local overabundance of massive progenitors and can remove the 
difficulty of the deficit of relatively bright 0.1 ROSAT PSPC cts s~^) cooling neutron stars previously reported 
from models where only neutron stars from the Galactic disk were accounted for. 
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1. Introduction 

Observations of isolated neutron stars (INSs) are impor- 
tant for gaining deeper insight on their structure and ther- 
mal evolution and ultimately might prove decisive in un- 
veiling the physical properties of matter at ultra- high den- 
sities. Up to a decade ago, the only known INSs were 
active radiopulsars, with the addition of the 7-ray pul- 
sar Geminga. Despite the fact that X-ray emission from 
some radio-pulsars was already detected by Einstein, it 
was in the '90s that ROSAT, thanks to its sensitivity in 
the 0.1-2 keV band, gave a clearer picture of the faint X- 
ray emission produced by the cooling surface of the closest 
INSs. ROSAT, supplemented by more recent observations 
by Chandra and XMM-Newton, revealed a variety of types 
of behavior in the X-ray emission from INSs and its rela- 
tion with radio activity. 

In particular, a substantial contribution of ROSAT to 
this field has been the discovery of a group of seven radio- 
quiet, close-by, thermally-emitting INSs, dubbed ROSAT 
INSs (RINSs) and sometimes referred to as the "mag- 
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posed, either as conventional middle-aged cooling NSs or 
as very old NSs a ccreting the interstellar gas (see again 
iTreves et ai]l2000|) . Although no compelling evidence has 
been brought in favor of either picture as yet, it is now 
generally believed that at least two of these objects (the 
brightest ones, RX J18 56-3754 and RX J0720.4-3125, e.g. 
iBraie fc Romani2002(l are indeed relatively young cooling 
NSs. 

Since the seven RINSs have remarkably similar ob- 
served properties, it is quite natural to assume that they 
belong to the same class, i.e. all of them are close-by, cool- 
ing NSs. This, however, poses a major problem. A useful 
and standard way to study a population of sources is to 
compute the Log N - Log S distribution starting from a 
model, and then comp are it with observations. This has 
been done fo r RINSs bv lNeuhauser fc TriimpeJ l)l999(l and 
IPopov et ai]|2000b) . The main conclusion of the latter in- 
vestigation is that the typical spatial density of radio pul- 
sars in our Galaxy is too low to explain RINSs. In other 
words, the assumption that RINSs and ordinary radio- 
pulsars derive from the same parent population underpre- 
dicts the number of observed "coolers" , i.e. NSs which are 
hot and close enough for their thermal emission to be de- 
tected in X-rays. 

An obvious solution is to invoke a local (both in space 
and time) overabundance of NSs with respect to those 
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originating in the Galactic disk and seen now as radiopul- 
sars. The main goal of this paper is to investigate the 
possible origin of these objects. Here we suggest that the 
likely birthplace for many of the IN Ss in the Solar proxim- 
ity is the Gould Belt fsee lPopov et al. 2002 for a prelimi- 
nary discussion) . The Gould Belt is a compound of young 
stellar associations extending in a ring tilted at about 20° 
from the Galactic plane, roughly centered at the position 
of the Sun and with a radius of ~ 1 kpc. Originally discov- 
ered in the middle of the 19th century, the Belt was stud- 
ied in detail by Benjamin Gould (see Stothers, fc, Frog el 
Il974l for a historical overview). Most of the known close- 
by young star associations (in Ophiuchus, Orion, Perseus, 
Scorpius etc.) belong to the Belt. Since the Gould Belt is a 
relatively young stellar system it contains a comparatively 
large fraction of massive stars. As the Solar proximity is 
embedded within the Belt, it is likely to have harbored 
quite a number of massive young stars which are the pro- 
genitors of NSs. 

In Sect. 13 we perform a population synthesis of close- 
by young INSs, assuming that they are born both in the 
Gould Belt and in the Galactic disk with an assigned 
mass spectrum. The dynamical evolution of the young 
NSs population in the Galactic potential is then followed 
together with the thermal evolution. This allows us to 
compute the present surface temperature and spatial dis- 
tributions of young, close-by NSs and finally their X-ray 
count rate, once a model for interstellar absorption has 
been prescribed. This is used to compute the Log N - 
Log S distribution which is then compared with observa- 
tions in Sect. El We refer to ROSAT PSPC in evaluating 
the X-ray count rate and in confronting our model with 
data. Even if Chandra and XMM-Newton provide much 
better observations of individual sources and small fields, 
ROSAT all-sky data are still the most complete by far in 
soft X-rays. Discussion and conclusions follow in Sect.0] 

2. The Model 

In this section we describe the method we used to com- 
pute the local distribution of young, cooling neutron stars. 
The main ingredients of our model are: the spatial distri- 
bution of NS progenitors, the NS formation rate, the NS 
cooling history, and the interstellar medium (ISM) distri- 
bution (the latter determines the interstellar absorption 
and hence affects the X-ray count rate). Even if we are 
concerned with a young population (r ^ 10® yrs), its dy- 
namical evolution can not be neglected if the number of 
sources in a limited volume (size ^ 1 kpc) has to be as- 
sessed. In fact, a NS with a typical velocity ~ 300 km s"^ 
travels a distance ~ 300 pc in its lifetime as a cooler. While 
not much on a Galactic scale, such a displacement is non- 
negligible in evaluating the cooler population in the Solar 
proximity. For this reason we account for the dynamical 
evolution of NSs in the Galactic potential in our model. 

Our calculation proceeds in three steps. First, a spa- 
tial distribution of the progenitors is selected and the 
ratio of NSs born in the Galactic plane to those origi- 



nating from the Gould Belt is fixed, together with the 
birth rates (Sect. I2.1|l . Second, the dynamical evolution 
is followed assuming that NSs at birth receive a kick ve- 
locity drawn from a prescribed velocity distribution (see 
Sect. I2.2|l . Finally, we derive the X-ray flux as a function 
of age and position on the basis of an updated set of cool- 
ing curves fSect. IT^ . and translate this in ROSAT count 
rate for an assumed model of the interstellar absorption 
(Sect. 123 • 

We varied the parameters of our model (NS forma- 
tion rate, NSs mass spectrum, time of calculation, time 
step, spatial distribution). In the following subsections we 
focus on those parameters corresponding to the results 
presented in Sect. 3. 

2.1. Initial spatial distribution 

In our picture NSs arc born in the Galactic disk and in the 
Gould Belt with a constant rate over the entire calculation 
(see Fig. n for a schematic view). This is a reasonable ap- 
proximation since we are dealing with young objects, less 
than a few Myrs old. The Gould Belt is modeled as a very 
thin disk of 500 pc in radius and inclination to the Galactic 
plane of 18° with its center situated at 100 pc from the Sun 
in the Galactic anticenter direction. The cen tral region of 
the Gould Belt is devoid of massive stars (see PoodcI '19971 
for de tailed description of the Gould Belt, and Torra ct alj 
20001 for a shorter one). To mimic this, we assume that no 
NS is born in the central region of the Belt up to a distance 
of 150 pc from its center. 

We assume that the birthrate per unit area is inde- 
pendent of position, i.e. it is constant both in the Belt 
and in the Galactic disk. The two (constant) values are 
different, the one referring to the Belt being larger. To 
derive the NS birthrate, we rely upon direct counting of 
massive stars which are doomed to end in a supernova 
even t. This gives 2.9x 10~^^ massive progenitors per yr per 
pc^ ijTammann et al.lll994l) which implies a NS birthrate 
of 30 NSs Myr~^ in 0.6 kpc around the Sun, the region 
where data are available. In relating the NS and supernova 
birthrates we assumed that core collapse produces a neu- 
tron star in 90% of the cases. The above value gives a more 
reliable estimate of the local NS birthrate, at present, as 
compared with radio-pulsar and average supernova statis- 
tics. 

To get the separate contributions to the NS birthrate 
from the Belt and from the disk for the region inside 0.6 
kpc we proceed as follows. In one Myr, out of a total of 
thirty, twenty NSs are born in the Gould Belt accord- 
ing to the esti mate of 20-27 supernova events given by 
iGrenierl l|200fi) . The remainder is uniformly distributed 
in the Galactic disk in a ring extending from 100 to 600 
pc, to account for the known fact that the closest Solar 
neighbourhood is underpopulated by massive stars (e.g. 
^afz-Aocllaniz 2001). This is in agreement with the re- 
sults of iTorra et all l|200nl) . who estimate that about 2/3 
of the massive stars inside 0.6 kpc belong to the Gould 
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Fig. 1. R — z projection of the initial spatial distribution 
of NSs. The heavy lines mark the Gould Belt and the 
closer (< 3 kpc) part of the Galactic disk where NSs are 
assumed to be born. Note the two regions devoid of NS 
formation, in the Galactic disk with a radius of 100 pc 
around the Sun, and in the Belt with a radius 150 pc 
around its center. The 0.6 kpc sphere centered at the Sun 
is also shown for clarity. 

Belt, and the rest to the Galactic disk. Direct counting 
of massive progenitors within 0.6 kpc from the Sun ac- 
counts for all NSs originating in the Belt. However, the 
disk clearly extends well beyond 0.6 kpc from the Sun. 
Although NSs born far away are very dim their contribu- 
tion may become significant at low fluxes 0.1 ctss""'^). 
For this reason we decided to include also NSs born in 
the Galactic disk from 0.6 up to 3 kpc from the Sun. The 
NS birthrate per unit area in this region is assumed to be 
again constant and coincident with that in the disk inside 

0. 6 kpc. With this choice we find an agreement with exist- 
ing data on the supernova rate i n a region of 1 kpc around 
the Sun l)Tammann et al.lll994 . 

In summary, we adopt a uniform birthrate of 20 NS 
Myr^i in the Gould Belt and 250 NSs Myr^^ in the 
Galactic plane up to a limiting distance of 3 kpc from the 
Sun. Our approach underestimates the number of distant 

1 kpc) newborn NSs in the direction of the Galactic 
center because the rate of supernova events is likely to 
increase there. However, interstellar absorption towards 
the Galactic center is very high in the Galactic plane (see 
Sect. 12.4(1 . and this hinders the detection of sources at low 
Galactic latitude. We note that the local overabundance 
of young NSs is probably due not only to the enhanced 
number of massive stars harbored in the Belt but reflects 
also the fact that the star formation history in the Belt 
favors an enhanced supernova rate in the present epoch. 

2.2. Dynamical evolution 

A detaile d description of the e volutionary code may be 
found in IPopov et alJ l)2000albl) . Typically we calculate 
about a thousand evolutionary tracks in each run (up to 
~ 10^ per production run) and then normalize our results 
to the actual number of NSs born in the volume (~ 1000) 
during the 4.25 Myrs time interval, as discussed above. An 
evolutionary track is calculated for specified initial posi- 
tion and velocity of a NS. All tracks are used over their 
whole duration and are applied to NSs of different masses, 

1. e. for each track we have several different cooling histo- 
ries, the number of which depends on the number of differ- 



ent masses used in the run. As in Popov et al. (2000a, b), 
each track actually represents a population of NSs of dif- 
ferent masses continuously born at the specified location 
with a prescribed spatial velocity distribution, during the 
entire time interval of the calculation. 

NSs are assumed to be born with a MaxwcUian kick 
velocity distribution with an average of ^ 225 km s^^. 
Varying this parameter does not change our results signifi- 
cantly because we are probing young objects. The Galactic 
potential is assumed to be axisymmetric, resulting from 
the sum of three co ntributions: the disk, bulge and halo 
(see |Paczvnskilll990|) . The Sun is in the Galactic plane at 
a distance 8.5 kpc from the Galactic Center. 

2.3. Cooling curves and flux calculation 

To calculate the cooling of NSs we use the resul ts of 
the St. Peters burg group (see Kaminker et al.ll2002l and 
the review by lYakovlev et alJ 1999() . Cooling curves are 
for NSs of masses from 1.1 M0 to 1.8 M© with a step 
of 0.1 M0 (see Fig. Curves take into account all im- 
portant processes of neutrino emission . The equation of 
state (E OS) used in [K aminke r et all (|2P0^ was intro- 
duced by IPrakash et al.. 1,1988)) . More precisely, it corre- 
sponds to Model I of Prakash et al. for symmetry en- 
ergy and compression modulus of saturated nuclear mat- 
ter K = 240 MeV. The maximum NS mass in this model 
is 1.977 Mq. Since we base our calculation on a definite 
EOS, the star radius R is known once the mass M is fixed. 
For I.IM© < M < I.8M0 it is 12.2 km < R < 13.2km, 
where R is circumferential radius. The EOS used here cor- 
responds to matter composed of neutrons, protons, and 
electrons (no muons, hyperons, and no exotic particles). 
Neutron superfluidity was ignored. 

We start to compute the emitted flux when a NS has an 
age of 10000 yrs. This choice is motivated by the fact that 
the Vela pulsar (the youngest close-by thermally emitting 
NS) is slightly older than that and also because the NS 
birth rate used in our calculation corresponds to an event 
every ~ 10 000 yrs in the ~ 1.7 kpc region around the 
Sun. For each NS the calculation is stopped when the 
temperature drops below 10^ K. This happens at an age 
of 4.25 Myrs for the lightest NSs {M = 1.1 Af©) or less 
for more massive stars. NSs that cold could have been de- 
tected by ROSAT PSPC within a distance of only 10 pc 
with a count rate ~ 5 x 10^'^ cts s^^. We assume that 
emission comes from the entire star surface. This appears 
reasonable in the light of the relatively low pulsed fractions 
(^ 1 5-20%) detected so far in the majority of isolated NSs 
see iHaberl fc Zavlinll2002() . We do not take into account 
any reprocessing of the blackbody surface emission in the 
NS atmosphere. 

NS s are expected to hav e a mass spectrum (see discus- 
sion m IWooslev et alJ2002|) . In the calculations we present 
here NSs are taken to have a flat mass spectrum between 
1.1 and 1.8Mq. The mass spectrum (and clearly the mass 
of each star) are assumed to be constant during the calcu- 
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Log Age (yr) 

Fig. 2. Cooling: c urves for different NS masses 
ijKaminker et al.ll2002|) . Curves from top to bottom cor- 
respond to masses 1.1-1.8 Mq (step O.IMq). Here, and 
everywhere in the text, the temperature is the red-shifted 
surface temperature (i.e. that observed at infinity). 



For th e ISM distribution we use the same prescrip- 
tion as in lPopov fc Prokhorovl l)l998|) . The Local Bubble 
is modeled as a sphere with a radius of 140 pc and ISM 
density of 0.1 cm~'^. Typical column densities for sources 
inside the calculated volume are in the range A'h ~ 10"'^^- 
10^^ cm^. After the column density is evaluated for a cur- 
rent NS position, we calculate the unabsorbed flux cor- 
responding to the temperature, radius of the NS and its 
distance from the Sun, and apply the standard procedure 
to derive the ROSAT PSPC count rate. Outside the Local 
Bubble in directions close to the Galactic plane, absorp- 
tion starts to play a crucial role. That is why regions closer 
to the Galactic center (< 7 kpc), where the NS formation 
rate should be higher than in the Galactic disk in the 
Solar vicinity, cannot add many sources to our sample. At 
iVjj = 3 X 10^^ cm~^ even very young and hot low-mass 
NSs with kT « 0.1 keV cannot have been detected at a 
distance 1 kpc in the ROSAT Bright Survey (RBS) at the 
threshold of ~ 0.2 cts s~^. 

Except for the Local Bubble (and even that in a quite 
simplified way, see ISfeir et al ] ll999l for a more complete 
description), our model does not take into account small- 
scale irregularities in the ISM distribution. They can be 
important if one makes an attempt to produce a realistic 
map of RINSs distribution on the celestial sphere, but in 
the case of the all-sky averaged Log N ~ Log S distribution 
our approximation is adequate. 



lations. With the above values (maximum age 4.25 Myrs 
and time step 10000 yrs) each evolutionary track repre- 
sents 425 X 8 NSs of eight different masses born with a 
time step 10^ yrs at the same place with the same initial 
velocity in a period 4.25-0.01 Myrs ago. 

As one can see from Fig. 12 the main contribution to X- 
ray bright sources comes from NSs with masses < 1.5 Mq. 
Young isolated NSs which are observed in several SN rem- 
nants as compact X-rav sources al so should be relatively 
low-mass objects (jKaminker et al ] |2002I Observations of 
binary radio pulsars sugge st that NSs masses are strongl v 
peaked around 1.35 Mq l|Thorsett fc Chakraba t7 19991. 
We repeated our calculations using a similar distribution 
and found that it produces similar results, although the 
number of observable isolated NSs increases by ~ 30 %. 
Calculations which take into account the realistic mass 
function of massive progenitors in the Solar vicinity will 
be the subject of a future paper. 

2.4. ISM and absorption 

Since young cooling NSs are expected to emit most of their 
luminosity at UV/soft X-ray energies 20 — 200 eV, 
corresponding to temperatures ~ 10^-10^ K), interstel- 
lar absorption plays a crucial role with respect to their 
observability. Any attempt to derive the number of ob- 
servable cooling isolated NSs using the unabsorbed flux 
would result in a substantial overestimate. 



3. Results 

As discussed in the previous section, the dynamical evo- 
lution of young neutron stars, together with the calcula- 
tion of their X-ray flux, allows us to derive the Log N - 
Log S for these sources. In particular, having fixed the 
number of NSs originating in the Gould Belt and in the 
Galactic disk, we can compute the separate contributions 
of these two sub-populations to the total Log N - Log S. 
Our main results are presented in Fig. |3| where the total 
(disk-|-Belt) and the disk-only distributions are compared 
with observations. All curves refer to the whole sky, i.e. the 
angular coverage used here is the entire solid angle. The 
observed Log N - Log S has been derived from ROSAT 
data of the seven RINSs and six other close-by young iso- 
lated NSs (see Table 1). Contrary to RINSs, the latter 
sources exhibit a composite spectrum with a non-thermal 
high-energy tail superimposed on the thermal component. 
The total count rate has been used in this case, but we 
stress that the non-thermal contribution is sizeable only 
for the Vela pulsar and PSR 1929-1-10. Poissonian errors 
are assumed in both observations and simulations. The 
computed curves were generated with a sufficiently large 
number of stars to reduce the statistical noise which is 
particularly severe at large flux (see i i2.2(l . Statistical er- 
rors are shown for the observational points only, being too 
small to be appreciated in the plot of the simulated dis- 
tributions. Uncertainties arising from the assumed model 
parameters are more difficult to assess and could introduce 
non- negligible deviations. Test runs with different choices 
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Fig. 3. All-sky Log N - Log S distribution: filled diamonds 
are the seven RINSs and open circles Geminga, the "three 
musketeers", PSR 1929+10 and 3EG J1835+5918. We 
also show the RBS limit (Schwope et al. 1999). Upper 
curve: NSs born in the Gould Belt and in the Galactic 
disk (total birth rate 270 NS Myr~^). Lower curve: disk 
only (birth rate 250 NS Myr^i). 

of key parameters show that differences are of the order 
of 50%. Symbols which show observational points (filled 
diamonds or open circles) correspond to the type of the 
faintest object (RINS or not) which contributes to the to- 
tal number at the specified count rate. 

Previous investigations have convincingly shown that 
cooling NSs with the same spatial density of ordi- 
nary radio-pulsars cannot explain the observe d Log N - 
Log S distribution for RINSs (Neuhiiuser fc Trtimperl 
I1999I rPopov et al.ll2000blp . In particular, it was stressed 
that an additional contribution to the local population 
of young cooling NSs should be invoked to explain the 
relatively large number of bright sources at fluxes <; 0.1 
cts s~^. The Gould Belt is a natural candidate to provide 
the missing NSs in the Solar proximity. The population 
synthesis calculations presented here strongly support this 
claim and account not only for the observed distribution 
of RINSs but also for that of other young, close-by isolated 
NSs observed by ROSAT. 

As can be seen from Fig.O the Gould Belt provides the 
major contribution to the local population of young cool- 
ing NSs and the theoretical prediction gives a very good 
fit to the data, within statistical errors. The contribution 
of NSs born in the Galactic disk is not very important 
at relatively large fluxes (it leaves rooms for only a few 
sources with count rate ^ 0.1 cts s~-^), but, as expected, 
becomes dominant at lower fluxes 0.01 cts s~^) where 
far away stars contribute most. Fig. 13 clearly shows that 
the curve referring to the Galactic disk alone is always be- 



low the observed points with high statistical significance 
and below the RBS limit. The worsening of the agreement 
at low fluxes (^ 0.05 cts s^^) is to be attributed to the 
incompleteness of the X-ray sample. The deficit of very 
bright objects (^ 1 cts s~^) may be attributed to chance 
statistics. 

Our calculations show that there can be about 10 — 50 
unidentified isolated NSs in the ROSAT All-Sky Survey 
(RASS) at a limiting flux of ^ 0.015 cts s"^ depend- 
ing on parameters of the model. The number of INSs 
in the RASS /BSC has been recently investigated by 
iRutledge et al.l l)2003l) . who found that at most 67 sources 
could have been detected at the 0.05 cts s~^ level and 
have escaped identification. Our results are well within 
this limit. Also there may be a few unidentified RINSs at 
fluxe s ^ 0.1 cts s~^ at low Galactic latitudes (see also 
Schwope et al.lll999tl . Most sources should be observed 
at ± 20° from the Galactic plane towards the directions 
of lower absorption. Some of them may be Geminga-like 
objects with counterparts among unidentified gar nma-ray 
sourc es (also connected with the Gould Belt, see iGrenien 
200(t . Identification of these objects may prove important 
for constraining cooling models and the NS mass spec- 
trum. 

Absorption, the flat geometry of NS initial distribu- 
tion and the finite extension of the Gould Belt naturally 
explain the very fiat (slope flatter than -1) Log N - Log S 
curves in Fig.|31 Our model contains some degrees of free- 
dom (e.g. the NS mass spectrum, details of the formation 
rate and surface emission, a simple blackbody was used 
here) that can be varied. We nevertheless believe that even 
this simple picture (albeit based on realistic assumptions) 
gives a quite satisfactory explanation of the observed prop- 
erties of isolated NSs in the Solar proximity. Preliminary 
calculations which take into account realistic mass spec- 
trum of NS progenitors, atmospheric effects and different 
spatial and velocity distributions are in progress and show 
no qualitative changes. 

4. Discussion and conclusions 

At present about 20 nearby (distance < 1 kpc), young 
(age < 4.25 Myrs) isolated NSs are known (see Table 
The local NS population includes objects with differ- 
ent properties: radio-quiet, thermally-emitting NSs (the 
seven RINSs), Geminga and the Geminga-like object 3EG 
J1835-f5918 (these are probably active pulsars with ra- 
diobeams missing the Earth), radio-pulsars with observed 
thermal X-ray emission (PSR 1929-1-10 and the "three 
musketeers": the Vela pulsar, PSR 0656-^14, PSR 1055- 
52) and seven other pulsars which are not detected in 
X-rays. The latter objects are relatively old in compar- 
ison with the others and lie further away, as can be seen 
from Table ^ where the wide gap in the estimated age 
of PSR 1055-52 and PSR J0056-H4756 is apparent. For 
ages of a few Myrs, even very low-mass (and hence slow- 
cooling) NSs are too cold by now to have been detected 
with ROSAT at distances ^ 0.5 kpc. The case of PSR 
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Table 1. Local (r < 1 kpc) population of young (age < 4.25 Myrs) isolated neutron stars 



Source name 


Period, 


Count rate, 


P 


Distance 


Age'' 


Ref. 




s 


ROSAT cts s"^ 


10-15 ss-i 


kpc 


Myrs 




RX J 185635-3754 




3.64 




0.117"* 


~ 0.5 


[1,2] 


RX J0720.4-3125 


8.37 


1.69 


~ 30 - 60 






[1,3] 


IRXS J130848.6+212708 (RBS 1223) 


10.3 


0.29 


< lO"*? 






[1,4] 


RX J1605.3+3249 (RBS 1556) 




0.88 




— 


— 


[1] 


RX J0806.4-4123 


11.37 


0.38 


— 


— 


— 


[1,5] 


RX J0420.0-5022 


22.7 


0.11 


— 


— 


— 


[1] 


IRXS J214303.7+065419 (RBS 1774) 




0.18 








[6] 


PSR B0633+17 (Geminga) 


0.237 


0.54'= 


10.97 


0.16'' 


0.34 


[7] 


RX J1836.2+5925 (3EG J1835+5918) 


— 


0.015 


— 


— 


— 


[8] 


PSR B0833-45 (Vela) 


0.089 


3.4^= 


124.88 


0.294"* 


0.01 


[7,9,10] 


PSR B0656+14 


0.385 


1.92^= 


55.01 


0.762'= 


0.11 


[7,10] 


PSR B1055-52 


0.197 


0.35" 


5.83 


~ 1" 


0.54 


[7,10] 


PSR B1929+10 


0.227 


0.012"= 


1.16 


0.33'* 


3.1 


[7,10] 


PSR J0056+4756 


0.472 




3.57 


0.998'= 


2.1 


[10] 


PSR J0454+5543 


0.341 




2.37 


0.793'= 


2.3 


[10] 


PSR J1918+1541 


0.371 




2.54 


0.684'= 


2.3 


[10] 


PSR J2048-1616 


1.962 




10.96 


0.639'= 


2.8 


[10] 


PSR J1848-1952 


4.308 




23.31 


0.956'= 


2.9 


[10] 


PSR J0837+0610 


1.274 




6.8 


0.722'= 


3.0 


[10] 


PSR J1908+0734 


0.212 




0.82 


0.584'= 


4.1 


[10] 



") Radio-pulsar ages are estimated as P/(2P), 

for RX J 1856-3754 the age estimate comes from kinematical considerations QVate^^JLattim^r ^00^ . 
*") Distance to PSR 1055-52 is uncertain (~ 0.9-1.5 kpc) 
"=) Total count rate (blackbody + non-thermal) 
'') Distances from parallactic measurements 
*=) Distances from the dispersion measure 

(1) Treves et al. (2000), (2) Kaolan, van Kerkwiik & Anderson (2002), (3) Zane et al. (2002) . 

(4) Hambar yan et al. (2001), Habcrl (2003), (5) Habcrl & Zavlin (2002), 

(6) Za mpieri et al. (20011 . (71 Becker fc Trumper (1997 ). (8) Mirabal & Halpcrn (2001), 

(9)^^^o^^^l^ (20^^, (10) ATNF Pulsar Catalogue, ( http: / /wwwatnf .atnf.csiro.au /research / pulsar/catalogue/t 



1929-1-10 is intermediate since it is ~ 3 Myrs old but rel- 
atively close (~ 300 pc) and its X-ray emission is mainly 
non-thermal. 

As has been discus s ed by iNeuhaiiser fc Triimperl 
figgg") and IPopov et all l)2000b(l . the number of X-ray 
bright INSs is too large to be explained in terms of the av- 
erage density of radio-pulsars in the Solar neighborhood. 
This motivated the suggestion that a sizeable fraction of 
NSs never become active radio emitters, as suggested at 
the same time by considerations on young NSs in super- 
nova remnants llCxotthelf fc VasishtL 1200(1 albeit the ori- 
gin of the local NS population was left open. In this paper 
we have shown that thermally-emitting INSs are naturally 
explained as cooling NSs born in the Gould Belt. The rela- 



tively high local spatial density of young NSs is then due to 
the large number of massive progenitors in the young stel- 
lar associations which constitute the Gould Belt. Our anal- 
ysis lends support to the idea that RINSs (the "magnifi- 
cent seven") represent the slowly coolin g NSs and hence 
a clea n sample of NSs with Al ^ 1.3M0 l)Kaminker et alJ 

The computed Log N -Log S distribution for X-ray 
thermally-emitting INSs born in the Gould Belt and (to a 
lesser extent) in the Galactic disk accounts for all bright 
INSs in the Solar vicinity (see Fig. ^ and leaves room 
for at most 1-2 undetected sources above ~ 1 cts s~^. 
The absence of sources brighter than RX J1856 (Log S> 
0.56) is consistent with our calculations at the 2 a level. 
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Having neglected very young NSs (age < 10000 yrs) does 
not change our conclusions significantly. Only young stars 
with M < 1.35 Af0 (about 1/2 of their total number) 
have temperatures about 2 x 10^ K (see Fig. 2), which 
corresponds to a luminosity ~ 10'^'* erg s~^. With our 
NS formation rate we expect to have one such NS within 
a distance of 2.5 kpc. About half of the NSs born in 
that region should originate inside 1.7 kpc. Placing the 
source at this distance its unabsorbed flux is ^ 3 x 10^^^ 
erg cm~^ s~^. With a typical column density of 10^^ cm~^, 
Web-PIMMS gives - 1.3 cts s'^ for ROSAT PSPC. We 
can then conclude that stars younger than 10000 yrs can 
contribute only to fluxes < 1.3 cts s'-^. This translates 
into a ~ 20% difference in the Log N - Log S curve shown 
in Fig. 121 in the worst case. Typical uncertainties of the 
present model are at about the same level. 

The local (^ 1 kpc) INS population was estimated to 
comprise total about one hundred objects younger than 
~4 Myr, taking into account that some NSs born in- 
side 1 kpc can leave the local space in their lifetime. 
These NSs are not detected as radio pulsars, but tens 
of them could be identified in ROSAT surveys as dim 
sources. The beaming effect can be responsible only for 
pa rt of these young N S s to be radiosilent. According 
to iTauris fc ManchesteJ l)l998|) in fact, about 50-70% of 
young pulsars are not v isible from the Earth (see also 
iBrazier fc Johnstonll 19991) . In our model this would imply 
the presence of at least ^ 30 active pulsars with an age 
<j 4 Myrs in the Solar vicinity. However, only about 1/3 of 
such pulsars are observed (see again Table 0. A bimodal 
velocity distribution, with the high charctcristic velocity 
^ 500 km s^i (Arzoumanian. Chernoff & Cordes 2002), 
would reduce the discrepancy since a number of young pul- 
sars leave the local volume considered here. However, even 
barring observational bias against low-luminosity radio- 
pulsars, our model seems to support the argument by 
ICotthelf fc Vasishd l|200nl ). that at least half of the ob- 
served young neutron stars follow an evolutionary path 
quite distinct from that of the Crab pulsar. 
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